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Abstract 

We construct a family-unified model on a Z2 x Z2 orbifold in five dimensions. The model is 
based on a supersymmetric SU(7) gauge theory. The gauge group is broken by orbifold boundary 
conditions to a product of grand unified SU(5) and SU{2) x U(l) flavor symmetry. The structure 
of Yukawa matrices is generated by an interplay between spontaneous breaking of flavor symmetry 
and geometric factors arising due to field localization in the extra dimension. 
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I. INTRODUCTION 



The success of gauge coupling unification in the Minimal Supersymmetric Standard Model 
(MSSM) suggests there exists a Grand Unified Theory (GUT). In GUTs, one generation of 
fermions can be incorporated in one or more representations of a simple GUT gauge group. 
However, GUTs do not explain why there are three different families and do not shed 
any light on the pattern of the observed fermion mass spectrum and mixing angles. An 
immediate idea is to also assign a symmetry group for the generations, the so-called flavor 
group. Theories incorporating this idea can be divided into two broad categories. First, 
theories in which the flavor group and the GUT group are orthogonal. Second, theories 
that unify the flavor group and the GUT group into a larger simple group. The second 
approach is known as family unification Within the first category, realistic theories with 
continuous flavor groups SU(3), SU{2) x £7(1), U{1) as well as several discrete groups have 
been considered. Family unified models push the unification idea a step further and are 
esthetically more attractive. However, four-dimensional models of family unification usually 
suffer from the problem of mirror families, see Ref. for a review. 

In the last few years the notion of symmetry breaking by orbifold boundary conditions in 
extra dimensions has been revitalized. Orbifold breaking has been used to address various 
problems ranging from electroweak symmetry breaking [jjj], supersymmetry breaking to 
GUT model building For example, in Ref. (if orbifold compactification is used to break 

the SU(5) GUT group to the Standard Model group and solve the doublet-triplet splitting 
problem that is difficult to overcome in 4D models. One of the reasons we use orbifold 
boundary conditions is to give large masses to mirror fermions, as noticed for example in 
Ref. 0. 

In this article, we construct a family-unified model in 5 dimensions. To maintain the 
unification of the gauge couplings our model incorporates supersymmetry. The fifth dimen- 
sion is compactified and we impose orbifold boundary conditions on all fields propagating in 
the fifth dimension. The role of the boundary conditions is threefold. We use the orbifold 
breaking to get rid of mirror families, break family-unified gauge group to a product of GUT 
and flavor symmetry, and also reduce the amount of supersymmetry to M = 1 in 4D. Given 
that the quarks of the third generation are a lot heavier than the quarks of the first two 
generations, it seems natural that the light families form a doublet, while the third family a 
singlet under the flavor group. Consequently, models using SU(2) x [/(l) flavor symmetry 
are quite successful in reproducing the mass spectrum [8|. We embed the flavor SU(2) x U (1) 
and SU(^) GUT group in an SU(7) family unified gauge group. A similar setup was studied 
in Ref. [9|, where an SU(7) family unified model was considered. However, in Ref. |9j the 
GUT group is flipped SU(5) and the emphasis is on the doublet-triplet splitting problem. 
A number of authors discussed flavor in extra dimensions, see Ref. [7[ and references within. 

The SU(7) gauge group is broken by the boundary conditions to SU (5) x SU (2) x 
£7(1). Both the SU(5) and the flavor groups are broken further by expectation values of 
Higgs fields. The pattern of Yukawa matrices is generated by both spontaneous breaking 
of the flavor group and geometric factors due to field localization. Some of the fields in 
our model propagate in the bulk, while others are localized at the orbifold fixed points. 
Therefore, wavefunction overlap suppresses certain couplings with respect to others. All 
of the flavor physics takes place at very high energy scales, comparable to the GUT scale. 
Supersymmetry breaking terms are of order the electroweak scale and are irrelevant for the 
discussion of flavor. We will not discuss the breaking of M = 1 supersymmetry in any 
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detail since this is not the focus of this paper. Any standard mechanism of communicating 
supersymmetry breaking in a flavor-diagonal manner could be incorporated into our model. 
Standard gauge mediation could operate if SUSY breaking and messenger fields are 
localized at an orbifold fixed point. By extending the model to one more dimension one could 
create an appropriate setup for either anomaly mediation 11 1 or gaugino mediation 0] of 
supersymmetry breaking. 

In the next section, we describe the field content and interactions needed to produce 
Yukawa matrices. We summarize our results in Section lTTTl The details concerning numerical 
determination of the high-energy parameters from the data are presented in Appendix [B] 



II. THE MODEL 

Our model is based on a supersymmetric field theory in five dimensions. The fifth dimen- 
sion is compactified on a (Z2 x Z2) orbifold. We parameterize the fifth dimension, described 
by coordinate y, as an interval with y G [0, ^]. This interval can be thought of as obtained 
from a circle [0, 27ri?] by identifying points related by reflections around two perpendicular 
axes. Under these reflections, y ~ — y and y ~ 71 — y such that the circle is equivalent to the 
y G [0, ^] interval. We denote these reflections as P and P', respectively. 

An arbitrary bulk field configuration can be decomposed into the eigenstates of the reflec- 
tions P and P' . Since P 2 = P' 2 = 1 the eigenvalues must be ±1. Of course, the eigenstates 
of the reflections have either the Dirichlet or Neumann boundary conditions at the end 
points of the interval. The Kaluza-Klein (KK) decomposition of a bulk field 0(2^,?/) into 
four dimensional mass eigenstates can be classified according to the two parities: 

y) = E - 7 =±=<f>^) cos 2 -f, (i) 

n=o V2-W.R ti 



00 



1 , 2 n+i/^„„„(2n+l)y 



x",y) = Y,-r^K- "V)cos 




R 

[x*,y) = £-^^+V)sii/"" 



R 



R ' 

where x M is the four dimensional coordinate and the subscripts refer to the parities under 
the P and P' reflections. The five dimensional Lagrangian has simple dependence on y when 
the fields are expressed in terms of KK states. The integral over the fifth dimension can 
be performed explicitly. One obtains then a four-dimensional Lagrangian describing a KK 
tower of four dimensional fields. The KK states specified in Eq. (0) have masses 2 "^j~ , 
and respectively. The only massless 4D field is 0° + (a; M ). 
It turns out that the compactification scale in our model will be comparable to the GUT 
scale. The massive states will therefore be too heavy to correspond to observable states. 
The fields of the MSSM will come from the zero modes of the KK decomposition, as well as 
from brane fields localized at the endpoints of the interval. 
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A. Fields and interactions 



We now begin to describe our model in detail. The 5D bulk theory is an M = 1 SUSY 
theory with an SU(7) gauge group. Such a theory has 8 supercharges and corresponds 
to M = 2 SUSY in four dimensions. However, the boundary conditions preserve only 4 
supercharges, so that below the compactification scale the theory is a four dimensional 
J\f = 1 theory. 

There is an arbitrary choice of how the reflection symmetry is represented in the 
space of gauge transformations. We choose the action of the two parities on the fun- 
damental representation of the SU(7) group to be P = diag{l, 1, 1, 1, 1, 1, 1} and P' = 
diag{ — l, —1, —1, —1, —1, 1, 1}. Consequently, an arbitrary tensor representation of SU(7), 
<f>ki"_, transforms as 

v) = v^p?pfp; l '---&:M, 

under the two parity transformations, where 77^,77^ = ±1 are the overall, "internal", parity 
eigenvalues. For a free field the parities can be chosen arbitrarily. Interaction terms correlate 
the parities of different fields. For example, the invariance of the supersymmetric Lagrangian 
imposes relations between parities of different components of superfields. 

The 5D gauge multiplet contains a vector Am, two gauginos Ai, A2, and a real scalar E, 
all of which transform in the adjoint representation of SU(7). We use the upper case Latin 
letters to denote 5D Lorentz indices, and the lower case Greek letters to denote 4D indices. 
The 5D SUSY Lagrangian is invariant under the reflections if 

V A , = ~Va 5 = ~Vv, V Xl = -V X2 , (2) 

as well as an identical set of relations for P' . We choose i] A = tj'a^ = V\ 1 = v' x — 1- 

Upon compactification, the first reflection breaks the M = 1 5D SUSY to M = 1 4D 
SUSY since both A 5 and A 2 obtain large masses. Meanwhile, A^ and Ai contain the zero 
modes that transform exactly as the 4D M = 1 vector multiplet. Since we do not embed the 
parity transformations into the R symmetry Af = 1 supersymmetry in 4D is preserved. The 
second reflection breaks the gauge group from SU (7) to its SU (5) x SU (2) x U (1) subgroup. 
More precisely, on the brane located at y — ^ the gauge group is broken, while in the bulk 
and on the brane y — the full symmetry remains. Besides the gauge multiplets, we put the 
5D hypermultiplets in the bulk. Under the SU (7) symmetry, the hypermultiplets transform 
as 1 + 7+35 + 21. 1 A hypermultiplet corresponds to two 4D chiral superfields with opposite 
parities ^ c }: 

r]q, = -r]q,c, (3) 
= -r^c (4) 

We choose 77^ = 77^ = 1 so that all the massless fields come from \1>. These massless fields can 
be expressed in terms of representations of the unbroken gauge groups SU (5) x SU (2) x U(l). 



1 The SU(7) field content coincides with an 5*0(14) spinor 64 when the spinor is written in the SU(7) basis. 
This suggests that our model may be embedded in a larger symmetry group. 
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These representations are T a (10, 2)_i, F a {5, 2) 3 , S ,a (l,2)_ 5 , and a neutral field (1, 1) , 
where a = 1, 2 is the SU(2) index. This set of 4D massless fields is free of gauge anomalies. 
In addition, there is no 5D anomaly either in the bulk or on the branes k| • We interpret the 
SU (5) as the GUT group and the SU{2) x U{1) as a flavor group. The zero modes (10, 2)_! 
and (5, 2) 3 are chosen to be the light two families of fermions with their superpartners and 
(1, 2)_ 5 might be the right handed neutrinos with their superpartners. 

On the asymmetric brane we add all other superfields that are necessary to complete 
the MSSM spectrum and break the GUT and flavor symmetries. Since the SU(7) gauge 
symmetry is broken on this brane, the localized fields do not need to form complete SU(7) 
multiplets. First, we choose the third family to be SU(2) singlets: T 3 (10, l)o, -^(5, l)o- 
Second, the SU(5) symmetry is broken by the conventional Higgses: £(24, 1) , H(5, 1) , 
H{5, l)o, K(45, l)o, and K(45, 1)q. We need to include the 45 and 45 representations to 



avoid the undesirable relations m 



1-1 



m. 



and m e ~ at the GUT scale. In addition, 



we introduce the following fields to break the flavor SU(2) x U(l) symmetry: 0T,a(l,2)i, 
T (l,2)_i, 0j? a (l,2)_3, F (1,2) 3 , ^(1, 1)_2, and ip(l, 1)2. We will refer to these fields as 
"flavons". We summarize the light field content in Table 1. 





bulk fileds 


brane fields 


matter 


r°(10,2)_i, F a (5,2) 3 , 5 a (l,2)_ 5 


T 3 (10,l) , F 3 (5,l) 


Higgs 
fields 




£(24, l)o, ff(5,l)o,i7(5,l) 
^(45,l)o, K(45,l) 


flavons 




0T,a(l,2)i, ^r(l,2)_i, V(l,l)-2,?(1,1)2 
^, (l,2)_ 3j 0^(1,2)3 



TABLE I: Light chiral superfields and their 577(5) x SU{2) x U(l) charges. 

Having listed the field content of our model, we now show how to obtain appropriate 
Yukawa couplings. Since our model is five- dimensional, the underlying field theory is non- 
renormalizable and has a cutoff A that is roughly two orders of magnitude larger than the 
compactification scale. To specify the Yukawa couplings we write the superpotential in terms 
of four-dimensional fields that is the brane fields and the zero modes of the bulk fields: 



W 



T 3 T 3 H + T 3 F 3 H + 



A 



T 3 T(j) T H + F 3 T<PtH + T 3 F<j) F H + T 3 T<f) T K 

1 



+F 3 T<p T K + T 3 F<p F K + TFipH + TF^K 
+T(f) T F(j) F H + T(f) T T(f) T K + T(j) T F(f) F K + TTZHip + TTZKip 



+ —[T<p T T<p T H 



(5) 



This superpotential is valid for a 4D theory below the compactification scale, so it is easy 
to keep track of dimensions of operators. We have explicitly indicated the 4- suppression of 
dimension five and dimension six terms. Different terms in Eq. (JSJ) are not related by any 
symmetries, so each term comes with a different coefficient. We have omitted the coefficients 
of operators for now. We will define and determine these coefficients in Appendix IB1 

In Eq. © we included almost all dimension five and six terms allowed by the gauge 
symmetries. We have omitted the couplings of K and K to the third family fields T 3 
and F 3 . Also, the flavons <fi F , 4> T do not appear in Eq. (JSJ. Since the superpotential is 
not renormalized it is technically natural to exclude certain terms. However, we can assign 
global symmetries to our fields such that the unwanted terms in Eq. (JSJ are prohibited. Such 
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symmetries also prevent (f) F and <f> T from appearing in the higher order terms, like i, -h, 
etc. For any term, the addition of the 577(5) adjoint £ is allowed by gauge symmetries. As 
we will show later, (S)/A is small. Thus, we also omitted terms with powers of £ whenever 
they would modify a lower order term that is more important. 

For the sake of clarity, we have also omitted geometric suppression factors in Eq. 
These factors are important for the structure of Yukawa matrices and are written explicitly 
in Sec. Ill Bl Such factors arise when bulk fields couple to brane fields because bulk fields 
propagate in 5D and their overlap with brane fields is small. A geometric factor is p = 
-j==. For a given term, the number of powers of p suppressing the term is equal to the 

number of bulk fields present in the term (lij . 

The 45 representation, K, and its conjugate K contain one SU (2) doublet each. Together 
with the doublets coming from H and H there would be four light doublets. We assume 
that one linear combination of doublets from H and K gets a large mass and the orthogonal 
linear combination remains light. The same thing takes place for H and K. We outline how 
to realize this in Appendix |XJ We denote the light mass eigenstates as h Ut d and the heavy 
ones as h' u d . In terms of these mass eigenstates 

H u = smah u + cos ah' u , K u = cos ah u — sin ah' u , (6) 
H d = sin jhd + cos jti d , K d = cos ^h d - sin jh' d , (7) 

where H u d and K u d represent the SU{2) doublet components of the corresponding fields. 
The mixing angles a and 7 are free parameters. For convenience, we define v = cot 7 and 
v' = cot a. 

We assume that the flavons and S get D-flat, SUSY-preserving, VEVs: 

<<M fa) fo\ (<M @ F ) fo\ (V) 65) 

a, (SJ 



A A \e J 1 A A \e' J' A A 

M = Diag{~ < y ) -|tf,-| < y, < y ) tf}. (9) 

Supersymmetry is only broken by weak-scale soft masses. We do not specify the superpo- 
tential that produces these VEVs as it is not essential for our discussion, but it would not 
be difficult to do so. 



B. Yukawa matrices 

Using Eqs. (0) through (0) it is straightforward to write the Yukawa matrices in terms of 
the MSSM superfields. These couplings arise at the GUT scale after (E) breaks 577(5) to the 
Standard Model. We denote the Yukawa matrices as Y u , Y d and Y\ for the up quarks, down 
quarks, and charged leptons, respectively. We do not consider the neutrino mass matrices 
or CP violating phases in this article. The Yukawa matrices are given by 



17 







p 2 aS(l 



-p 2 a5(l + v) p 2 e 2 (l + v) 
\pt{l+v) 





\pe(l + v) 
1 



sin 7, 



(10) 



6 



Ya 



Y, 



V 



p 2 a(l + v') \ 

-p 2 a(l + v') p 2 ee'(l + v') pe(l + v') sin a, 
pe'(l + <) 1 j 

p 2 a(l - 3v') \ 

-p 2 (7(l-3^/) p 2 ee'(l - 3i/) pe(l - 3v') sin a. 
pe'(l-3t/) 1 j 



(11) 



(12) 



We indicated the matrices with the proportionality sign because in the interest of clarity 
we omitted arbitrary coefficients that are also missing in Eq. (J5J). Compare Eqs. (|B2|) - (|B4|) 
in Appendix [B] that contain the full set of coefficients. The factors of | that appear in the 
(2,3) and (3,2) elements of Y u arise because the terms T 3 T<PtH and T^T^K contribute to 
both elements. 

We determine the magnitudes of the elements of the Yukawa matrices at the GUT scale 
by using the renormalization group equations for these matrices and comparing them with 
the masses and the CKM angles at the weak scale. The fitting procedure is described in 
Appendix |Bj As we will see the experimental data can be fitted quite accurately. Before we 
present the results let us make several comments about the matrices (|10 p -([12j ) . 

If e ~ e' the structure of the quark Yukawa couplings is very similar to the 4 texture zero 
symmetric quark mass matrices discussed in the literature. See Ref. 0] for a review. (Since 
the matrices discussed in Ref. are symmetric the off-diagonal zeros are the same element 
and counted as one zero.) The only difference is that our matrices are antisymmetric in the 
(1,2) and (2,1) indices 0, [B . The matrix for the leptons is similar to that for the down 
quarks. The zeros in the Yukawa matrices ()10|) - (jl2}) are exact provided that <p F and <p T are 
absent in the superpotential in Eq. (J5J). 

As discussed in Ref. this kind of matrices can give us the approximate relations 



IK, 
Vu 

V ts 







V ™ c 


V m s 






V m s' 





(13) 
(14) 



where is a CP violating phase defined in Ref. [17j. Of course, our matrices have only 
real elements, so the phase in the first relation is absent. In order to avoid the undesirable 
relation 



ub 



Vr 



cb 




(15) 



the (2,3) and (3,2) elements of the Yukawa matrices should be sufficiently large-much larger 
than m s /mb and m c /m t for d- Yukawa matrix and u- Yukawa matrix, respectively. 

The (2,3) and (3,2) elements of Y d and Y\ are too large to be neglected compared to the 
(2,2) and (3,3) diagonal elements. The off-diagonal elements are different for Yd and Yf. they 
are pe(l + v') and pe(l — 3v'), respectively. Thus, the off-diagonal elements affect the largest 
eigenvalue of the matrix differently for the bottom Yukawa and the r Yukawa. Therefore, the 
bottom-r unification is not exact and the b and r masses can be fitted accurately. Similar 
observation was made in Ref. 

It is not possible to uniquely determine all the parameters in Eqs. firU|) - (fT2"j) like e, v, 
etc. because there are arbitrary coefficients in front of every term, see Appendix El 
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Only certain combinations of the coefficients at and other parameters appear in the Yukawa 
matrices. We would like all coefficients to be close to one since they are dimensionless 
couplings. We get the values of a, to be close to one by choosing the remaining parameters 
as follows 

pe = ^ pe ' = w A = 2cM' 6 = ^ v ' = h v = l (16) 

We assume that p ~ 0.1 and infer the VEVs of flavons: e ~ 0.33, e' ~ 0.75, and a « 
0.15. These VEVs are smaller than the cutoff A, but e' is quite close to 1. Together with 
5 « 0.05 and the GUT scale M GUT ~ 2.8 x 10 16 GeV, we get A ~ 5.6 x 10 17 GeV and 

~ 8.8 x 10 15 GeV. These numbers give a resonable separation of the cutoff, the GUT, and 
the compactification scales However, the VEVs of the flavons are sufficiently close to the 
cutoff scale that higher dimensional operators may play an important role in the generation 
of Yukawa couplings in our model. The values of parameters in Eq. ()16|) correspond to 
tan (3 = 47, but we could make similar choices for other values of tan (5. 

What is interesting is that the flavons' VEVs: e, e', a are of the same order. This is 
very different from many 4D models where the flavons usually obtain hierarchical VEVs in 
order to produce hierarchy in the Yukawa matrices. The geometric suppression factor does 
contribute in our model to generating small ratios. 



III. SUMMARY 

The underlying theory for our model is a 5D SUSY theory with an SU (7) gauge group. 
Compactification of the fifth dimension on a Z 2 x Z 2 orbifold breaks SUSY to M = 1 in 4D 
as well as breaks SU(7) to GUT SU(5) times flavor SU(2) x U(l). The compactification 
scale is very close to the GUT scale, it is just a factor of three smaller than the GUT scale. 
Thus, our model is an ordinary SUSY GUT almost all the way to the GUT scale. In addition 
to symmetry breaking by boundary conditions we introduce two types of Higgs fields. First, 
standard Higgs fields that break GUT symmetry down to the Standard Model and give 
masses to the quarks and leptons. Second, flavon Higgs fields whose role is to completely 
break the flavor symmetry. The flavor symmetry is broken close to the GUT scale. 

Bulk multiplets contain zero modes corresponding to the two lightest families that trans- 
form as a doublet under flavor SU(2). The third family is a singlet under the flavor symmetry 
and it is localized at one of the orbifold fixed points. The SU(7) gauge symmetry is not 
preserved at the fixed point where the third family is localized. Therefore, the third family 
does not come from a complete SU(7) multiplet and is a flavor singlet. As far as the fla- 
vor symmetry and the light fields are concerned our model is very similar to the 4D model 
described in Ref. 

Our main goal was constructing a realistic pattern of Yukawa matrices at the GUT scale. 
We were only concerned with the quark and charged lepton sectors and completely neglected 
the neutrino sector. The Yukawa couplings come from the superpotential in Eq. (J5J), which 
we chose to resemble the "four zeros" texture described in Ref. [15] . The resulting Yukawa 
matrices, omitting a number of dimensionless constants of order one, are given in Eqs. ()10jl. 
(fTTfl. and flUD- 

The orders of magnitude of different elements of the Yukawa matrices are governed by 
three different effects. The first effect is the geometry of our model. The couplings that 
involve both localized fields and bulk fields are suppressed due to small wavefunction overlap 
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between such fields. Second, the SU(2) x U(l) flavor symmetry is broken by three different 
flavons and their conjugates. Among the three flavons there are two SU(2) doublets and one 
singlets. All flavons are charged under the U{1). We do not count separately the conjugates 
of the flavons because the VEVs of flavons with the conjugate quantum numbers are identical 
to maintain SUSY above the weak scale. Third, the up and down sectors are distinguished 
by the mixing of the Higgs doublets that come both from the 5 and the 45. The light up 
and down Higgs doublets come from different linear combinations of 5 and 45. Of course, 
any value of tan (3 other than 1 also differentiates the up and down sectors. 

Our model has too many free parameters to be predictive. What we accomplished, 
however, is generating the Yukawa matrices in terms of a few small parameters: flavon 
VEVs, defined in Eqs. (JBJ) and (|16p. and the geometric suppression factor. By matching 
to the observed fermion mass spectrum and quark mixing angles we determined the 13 
nonzero parameters in the Yuklawa matrices, see Eqs. (jB5|) - (jB7|) . We chose the undetermined 
parameters such that the dimensionless couplings are close to one. 

What is interesting is that given a few arbitrary choices all dimensionless coefficients are 
of order one. Moreover, many of the coefficients listed in Eq. (jBll|) are very close to one. 
All the large ratios are determined in terms of the geometric suppression factor and a few 
flavon VEVs that are of the same order of magnitude. Obviously, a more fundamental and 
predictive structure of flavor is still missing. However, it is conceivable that the flavor could 
be generated from an interplay between geometry and flavor symmetries. 
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APPENDIX A: MIXING OF HIGGS DOUBLETS 

We briefly comment on the mixing of the Higgs doublets coming from the 45 and 5 
representations and their conjugates. As we explained in Sec. |Hj we need the 45 and 45 
representations to avoid the equality of the lepton-down quark Yukawa couplings in the two 
light families. The problem is similar to 50(10) unification, where one needs to introduce 
larger Higgs representation in addition to the 10-dimensional Higgs to incorporate realistic 
Yukawa couplings. The additional Higgs fields, for example 126, would produce too many 
light doublets. A simple solution was presented in Ref. 20]. Similar solution works in the 
SU(5) case and we outline it here for completeness. 

We supplement the Higgs fields H, H, K, and K introduced already by another pair of 
45 and 45 Higgs fields. Let us refer to the new fields as K x and K\. We assume that the 
superpotential for these Higgs fields is given by 

W m g gs = I^HH + HYJI + Hi:K 1 + ISiTj + M X K\K + M 2 KK 1 . (Al) 

In the equation above E is the £77(5) adjoint field that develops an SU(5) breaking VEV 
given by Eq. Q and (E) oc 5 A. Mi and M 2 are arbitrary mass parameters that are com- 
parable to the GUT scale. We also assume that « — 5A so that the SU(2) doublets in H 
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and H are light. The mass matrix for the doublets arising from H, H, K, K, K\, and K\ 
has the following structure 





(H d \ 


T 


( o 


5A 





\ 


/ H u \ 




M = 


K ld 




5A 





M 2 






(A2) 




\K d ) 




{ 


Mi 





/ 


V K u J 





where we assumed that p + 5A = 0. The light eigenvalues of this mass matrix are two 
doublets 

h u = cos a H u + sin a K u , (A3) 
h d = cosjH d + siwyK d , (A4) 

where tana = <5A/M 2 and tan 7 = 5 A/ Mi. There is no reason to assume that Mi and M 2 
are equal, so the mixing angles of the up and down Higgs doublets are, in general, different. 
Clearly, the remaining Higgs doublets have masses of order the unification scale and so do 
other components of K, K, K±, and K\. 



APPENDIX B: FITTING TO THE DATA 

We now describe our procedure for determining the Yukawa matrices at the GUT scale. 
As we mentioned in Sec. Ill A[ there are no relations between the various terms in Eq. (jSJ 
since different terms are not related by symmetries to one another. Therefore, one needs to 
include arbitrary coefficients of order one in front of every term. After including the missing 
coefficients Eq. (jSJ) becomes: 



W = ai T 3 T 3 H + a 2 T 3 F 3 H + - a 3 T 3 T<j) T H + a 4 F 3 T<j) T H + a 5 T 3 F<f) F H + a 6 T 3 T<f) T K 



+ [a n T(f) T T(j) T H 



+a 7 F 3 T(j) T K + a 8 T 3 F<j) F K + a 9 TFipH + ai TFipK 
+a l2 T(f) T F(f) F H + a 13 T(j) T T(j) T K + auTfaFfoK + a 15 TTZHijj + a^TTZKip 
The corresponding Yukawa matrices are then 



Bl) 



Y,. 



Y, 



Y 







p 2 aS(a 15 + a m v) 







-p 2 a5(a 15 + a w v) p 2 e 2 (a 11 + a 13 v) f{a 3 + a 6 v) 
\ f(a 3 + a e v) ai J 

I p 2 cr(a 9 + a w v') N 

—p 2 a(a 9 + a w v') p 2 ee'(a 12 + a^v') pe(a 4 + a 7 v' 
V pe'(a 5 + a 8 v') a 2 J 

( p 2 cr{a 9 - 3a 1Q v') 

— p 2 cr(a 9 — 3a w v') p 2 ee'(ai 2 — 3a u v') pe(a 4 — 3a 7 v') 
\ pe'(o5 - 308^') a 2 



sin 7, 



sin a, 



vB2) 



sin a. 



(B3) 



(B4) 



The Yukawa matrices are defined in terms of 16 coefficients a i; 3 flavon VEVs, S VEV, 
geometric factor p, and two Higgs mixing angles: a total of 23 parameters. However, several 
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of our parameters only appear in particular combinations, which allows us to eliminate the 
" unobservable" combinations: 



/ ci \ 



u 



-Cl C 2 C 3 T], 



(B5) 



V c 3 1 J 
( c 4 \ 



-c 4 c 5 c 6 C> 



(B6) 



V c 7 1 J 
I c 8 \ 



-C 8 Cg Cio C- 

V en 1 j 



(B7) 



We are left with 13 parameters: C\ through en, 77, and £. The experimental data, gives nine 
masses and three real angles in the CKM matrix. Including complex phases in our Yukawa 
matrices would introduce too many free parametres, so we omit the phases. If we performed 
a fit with the phases present the values of the real parameters might change slightly, but 
such change would not affect the structure of the Yukawa matrices. 

To obtain the GUT scale values we use the following fermion parameters and the gauge 
couplings at the scale Mz jljj as inputs: 

at x = 0.016829, a 2 = 0.033493, a 3 = 0.118, 

m u = 2.33 ± 0.435 MeV, m c = 0.677 ± 0.0585 GeV, m t = 181 ± 13 GeV, 

m d = 4.36 ± 1.13 MeV, m s = 72 ± 23 MeV, m b = 3.00 ± 0.1 GeV, 

m e = 486.84727±0.00014keV,m /1 = 102. 75138±0. 00033 MeV, m T = 1.74669±0.000285 GeV, 

V us = 0.2205 ± 0.0018, V cb = 0.0373 ± 0.0018, \V ub /V cb \ = 0.08 ± 0.02. 

Given the structure of Yukawa matrices described in Eqs. (jB5|) - (|B7l) at the GUT scale, we 
use the one loop renormalization group equations in the MSSM |22j to compare with the 
weak scale data. We set the GUT scale to be Mqut = 2.80 x 10 16 where the three gauge 
coupling constants unify. The one-loop running of the gauge couplings does not involve 
the Yukawa couplings, so the gauge couplings are determined at all scales before fitting the 
Yukawa matrices. We neglect the fact that the first few KK modes appear below Mqut 
since RMqut ~ 3.2. Since the logarithm of RMqut is small we can neglect the effects of 
the KK modes below Mqut and we use 4D RGE equations. 

In practice, we numerically evaluate the RGE equations from the GUT scale down to the 
weak scale. We then compare the results of the RGE running with the data and evaluate the 
X 2 using the experimental errors. The errors are severely underestimated this way because 
threshold corrections and two-loop effects are much larger than the experimental uncertainty 
of the lepton masses. However, since we have more parameters than the number of inputs 
we are able to get a good fit. For example, we present the numerical fit for tan j3 = 47 below. 



u 



I -0.0001050 \ 

0.0001050 0.005335 0.05848 
v 0.05848 1 j 



1.053, 



(B8) 
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Y, 



Y 



0.004744 \ 
-0.004744 0.006898 0.1009 0.1597. 
0.2205 1 







-0.003774 







0.003774 -0.007916 -0.1526 
-0.3161 1 



0.4597. 



(B9) 



(BIO) 



The total \ 2 for this fit is 4.654. The \ 2 is dominated by the errors from rrid and m T , but 
none of the two masses with the poorest fit deviates by more than 1.5cx from the experimental 
value. 

To extract the physical parameters p, a and so on, we need to remember that the co- 
efficients ai,...,a 16 are close to 1. A choice of parameters is given in Eq. (jTUjl and the 
corresponding coefficients are 



ai = 1.266, 
a 7 = 1.019, 



a 2 = 0.8935, a 3 = 3, a 4 = 1.006, a 5 = 1.029, a 6 = 2.162, 
8 = 0.9590, a 9 = 1.557, aio = 0.7611, a u = 3, ou = 1.141, 



a i3 = 4.618, a u = 0.7942, 



ai5 



1, ai6 = 1.160. 



(Bll) 



Note that the coefficients 03, an and 015 in Y u are set by hand. There are too many free 
parameters to be uniquely determined from the 13 parameters in Eqs. (jB5|) - (JB7|) . so we need 
to arbitrarily choose some of them. 

Most coefficients eij are very close to 1 and certainly none of the coefficients deviate from 
one by an order of magnitude. There are many small additional contributions to this result 
that we neglected, for example higher-dimensional operators, threshold corrections, higher- 
loop effects. Since the coefficients tij are so close to multiples of 1 it is possible that the 
deviations could be accounted for by higher-order effects we neglected. This suggests that 
there could be a simple set of hidden symmetries responsible for this result. It certainly 
would be an exciting possibility. 
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